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We review the available experimental information on the thermodynamic properties of supercooled ordinary and heavy 
water and demonstrate the possibility of modeling these thermodynamic properties on a theoretical basis. We show that 
by assuming the existence of a liquid-liquid critical point in supercooled water, the theory of critical phenomena can 
give an accurate account of the experimental thermodynamic-property data up to a pressure of 150 MPa. In addition, we 
show that a phenomenological extension of the theoretical model can account for all currently available experimental 
data in the supercooled region, up to 400 MPa. The stability limit of the liquid state and possible coupling between 
crystallization and liquid-liquid separation are discussed. It is concluded that critical-point thermodynamics describes 
the available thermodynamic data for supercooled water within experimental accuracy, thus establishing a benchmark 
for further developments in this area. 



I. INTRODUCTION 

The peculiar thermodynamic behavior of supercooled wa- 
ter continues to receive considerable attention, despite several 
decades of work in this field. Upon supercooling, water ex- 
hibits an anomalous increase of its isobaric heat capacity and 
its isothermal compressibility, and an anomalous decrease of 
its expansivity coefficient.^ One explanation of these anoma- 
lies, originally proposed by Poole et al.,^ is based on the pre- 
sumed existence of a liquid-liquid critical point (LLCP) in 
water in the deeply supercooled region. The hypothesis of the 
existence of a critical point in metastable water has been con- 
sidered by many authors, as recently reviewed by Bertrand 
and Anisimov.-' In particular, several authors have made at- 
tempts to develop a thermodynamic model for the thermo- 
dynamic properties of supercooled water based on the LLCP 
scenario. ^"^ The existence of a liquid-liquid critical point in 
supercooled water is still being debated, especially in view 
of some recent simulations.^'^" The purpose of this article 
is to demonstrate that a theoretical model based on the pre- 
sumed existence of a second critical point in water is capable 
of representing, accurately and consistendy, all available ex- 
perimental thermodynamic property data for supercooled or- 
dinary and heavy water. While being conceptually close to the 
previous works by Fuentevilla and Anisimov^ and Bertrand 
and Anisimov,^ this work is the first comprehensive correla- 
tion of thermodynamic properties of supercooled water, in- 
corporating non-critical backgrounds in a thermodynamicaUy 
consistent way. 

This article is organized as follows. In Sec. 11 we review 
the currently available experimental information for the ther- 
modynamic properties of supercooled water. We also provide 
an assessment of the IAPWS-95 formulation' ' '^ for the ther- 
modynamic properties of H2O (designed for water at temper- 
atures above the melting temperature) when extrapolated into 
the supercooled region, confirming the practical need for an 
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improved equation of state for supercooled water. In Sec. Ill 
we formulate a thermodynamic model for supercooled water 
by adopting suitable physical scaling fields with respect to the 
location of a liquid-liquid critical point in supercooled water. 
In Sec. IV we show that this theoretical model yields an ac- 
curate representation of the thermodynamic property data of 
both supercooled H2O and D2O up to pressures of 150 MPa. 
In Sec. V we present a less-restricted phenomenological ex- 
tension of the theoretical model and show that this extension 
allows a representation of all currently available experimental 
data for supercooled H2O up to the pressure of 400 MPa. The 
article concludes with a discussion of the results and of some 
unresolved theoretical issues in Sec. VI and VII. 



II. REVIEW OF EXPERIMENTAL DATA 

Experimental data on the properties of supercooled water 
have been reviewed by Angell in 1982,'^'''* by Sato et al. 
in 1991,^^ and by Debenedetti in 2003.' Therefore, in this 
review, we focus on data published after 2003, and restrict 
ourselves to bulk thermodynamic properties. For a complete 
overview of the older data the reader is referred to the earlier 
reviews. Besides reviewing experiments, we also assess the 
performance in the supercooled region of the current reference 
correlation for the properties of water and steam, the "lAPWS 
Formulation 1995 for the Thermodynamic Properties of Or- 
dinary Water Substance for General and Scientific Use," or 
IAPWS-95 for short."-'^ Such an assessment has been car- 
ried out before,''^"'^ but not with all property data that are now 
available. Most of the data discussed here are provided in nu- 
merical form in the supplemental material.'^ 

A. Density 

Since the review by Debenedetti,' new data for the den- 
sity of supercooled water have been reported.^""^^ Most no- 
table is the recent work of Mishima,^' who measured the den- 
sity and compressibility down to 200 K and up to 400 MPa; 
see Fig. 1(a). More accurate density measurements have been 
pubUshed by Sotani et al.}'^ Asada et al.^^ and Guignon et 
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FIG. 1. (a) Location of the experimental H2O density data.'^"^^ The solid curve is the ice-liquid phase boundary;^^ the dashed curve is the 
homogeneous ice nucleation hmit.-^'*'-^^ The location of the dashed curve above 300 MPa is uncertain. At 0. 1 MPa and in the stable-liquid region, 
data from several older sources have been omitted for clarity, (b) Location of the experimental H2O density-derivative data. Ter Minassian et 
al?^ have measured the expansivity coefficient, other authors^ have measured the isothermal compressibility. 



ah, but their lowest temperature is 253 K, so in a larger tem- 
perature range Mishima's data are the only data available. The 
data of Guignon et al. are in the same range as the data of 
Sotani et al. and Asada et al. ; the maximum density difference 
between the data is 0.25%, which is within the experimental 
uncertainty. 

At atmospheric pressure, the density of supercooled water 
has been measured by several experimentalists.' We consider 
the measurements of Hare and Sorensen'** of 1987 the most 
accurate. They showed that their measurements were not af- 
fected by the 'excess density' effect, which occurs in thin 
capillary tubes and caused too large densities in their 1986 
experiments'^^ and in experiments of others. 

A comparison of the densities calculated from the lAPWS- 
95 formulation with the experimental density data of Hare and 
Sorensen,'^ of Sotani et al.}'^ and of Mishima^' is shown 
in Fig. 2. While the IAPWS-95 formulation reproduces the 
experimental density data at ambient pressure, the deviations 
from the formulation become larger and larger with increasing 
pressure. Especially at higher pressures, there is a sizable dis- 
crepancy between the IAPWS-95 formulation and the experi- 
mental data; the slope (or the expansivity) even has a different 
sign. 



B. Derivatives of tiie density 

Both the isothermal compressibihty and the expansivity 



coefficient of supercooled water have been measured; see 
Fig. 1(b). The most accurate compressibility data are from 
Kanno and Angell,'^^ whereas Mishima's^' data cover the 
largest temperature range. The only expansivity measure- 
ments are from Ter Minassian et al.^^'' Hare and Sorensen'**-^^ 
also have published expansivities (at 0. 1 MPa), but these were 
obtained from the derivative of a fit to their density data. 

According to Wagner and PruB,'^ the behavior of the ex- 
pansivity coefficient calculated from the IAPWS-95 formula- 
tion should be reasonable in the liquid region at low tempera- 
ture. However, from Fig. 3 we see that the IAPWS-95 expan- 
sivity is in error by up to 50% in the low-temperature region, 
even at temperatures above the melting temperature where the 
IAPWS-95 formulation should be vahd. The isothermal com- 
pressibihty calculated from the IAPWS-95 formulation agrees 
with the experimental data down to about 250 K and up to 
400 MPa, as shown in Fig. 4. However, at lower tempera- 
tures, the IAPWS-95 compressibihties do not even quahta- 
tively agree with the data. 



C. Heat capacity 

The isobaric heat capacity Cp of supercooled water has 
been measured only at atmospheric pressure.' Old measure- 
ments by Anisimov et al.^^ down to 266 K already showed 
an anomalous increase of the heat capacity at moderate super- 
cooling. In breakthrough experiments, Angell et a/.^'~^^ ex- 
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FIG. 2. Densities according to IAPWS-95 (curves). IAPWS-95 is 
valid to the right of the melting curve^-' Tm; the IAPWS-95 values left 
of the melting curve are extrapolations. The symbols represent exper- 
imental data of Mishima,^' Sotani et at. and Hare and Sorensen.'^ 
Mishima's data have been adjusted as described in Appendix C. 
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FIG. 3. Expansivity coefficient according to IAPWS-95 (solid 
curves: within region of validity, dashed curves: extrapolation). Sym- 
bols represent experimental data of Ter Minassian et al?^ and Hare 
and Sorensen.^^'-^' 



tended the range of measurements down to 236 K, and demon- 
strated that Cp keeps increasing with decreasing temperature 
(Fig. 5). More recent measurements by Archer and Carter,^"* 
also down to 236 K, do not perfectly agree with those of 
Angell et al?^ Archer and Carter suggest that the tempera- 
ture caUbration procedure of Angell et al. might cause a sys- 
tematic error in their measurements. Furthermore, Archer and 
Carter suspect that measurements of Tombari et al?^ (down 
to 245 K) could be affected by even more significant system- 
atic calibration errors. The measurements of Bertolini et al?^ 
down to 247 K agree with those of Angell et al. (after correc- 
tion, as described in Appendix C). 

Because there are no high-pressure measurements of Cp in 
the supercooled region, we mention here some experiments 
at high pressure P in the stable region. Sirota et aO'' mea- 
sured Cp up to 100 MPa and down to 272 K. Recently, Manya 
et al.^^ have measured Cp at 4 MPa from 298 K to 465 K. 
It turns out that the results of Manya et al. imply that the 
derivative {dCp/dP)T at constant temperature T is positive 
for pressures lower than 4 MPa, which contradicts the ther- 
modynamic relation {dCp/dP)r = -T{d^V /dT^)p, with V 



being the molar volume. Hence, the data of Manya et al. will 
not be considered in this paper. 

To our knowledge, the isochoric heat capacity Cy of super- 
cooled water has not been measured. The values presented by 
several investigators^^'^^'^^ were calculated from other ther- 
modynamic properties. 

The isobaric and isochoric heat capacities predicted by 
IAPWS-95 at 0.1 MPa are shown in Fig. 5. The Cp curve com- 
puted from IAPWS-95 follows the data of Angell et al.P to 
which it was fitted. Consequently, IAPWS-95 systematically 
deviates from the data of Archer and Carter,^"* as discussed 
by Wagner and PruB.'^ The isochoric heat capacity Cy does 
not exhibit a strongly anomalous behavior according to the 
IAPWS-95 prediction. 



D. Surface tension 

As listed by Debenedetti,' the surface tension of super- 
cooled water with vapor and with air has been measured by 
Hacker,'*" Floriano and Angell,^^ and Trinh and Ohsaka,^^ all 
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FIG. 4. Isothermal compressibility according to IAPWS-95 (curves). 
lAPWS-95 is valid to the right of the melting curve Tm; the IAPWS- 
95 values left of the melting curve are extrapolations. Symbols repre- 
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FIG. 5. Isobaric and isochoric heat capacity at 0.1 MPa according to 
IAPWS-95 (curves). Symbols represent experimental data.^^"^^ 



FIG. 6. Surface tension according to the lAPWS equation** 
(solid curve: within region of validity, dashed curve: extrapolation). 
Symbols represent experimental data from Hacker,^ Floriano and 
Angell,"*! and Trinh and Ohsaka.''^ 



face tension of liquid water in equilibrium with water vapor.^ 
For a comparison with experimental data, it is necessary to 
consider the conditions under which the surface tension was 
measured. Hacker,*" Floriano and Angell,*^ and Trinh and 
Ohsaka*^ measured the surface tension at atmospheric pres- 
sure in air. The influence of air at atmospheric pressure on 
the surface tension of water is usually neglected,*^ but for su- 
percooled water the effect may be significant since the sol- 
ubility of nitrogen and oxygen in water increases with de- 
creasing temperature. From high-pressure data*^ it is known 
that the surface tension of water with nitrogen or oxygen is 
lower than the pure-water surface tension. In Fig. 6 we show 
a comparison of the experimental surface-tension data with 
the values calculated from the lAPWS equation. As noted 
above, the experimental values of the surface tension suggest 
an inflection point at about 268 K. The extrapolation of the 
lAPWS equation does not show an inflection, and the differ- 
ence with Hacker's accurate data*" increases with decreasing 
temperature. The data of Floriano and Angell*' show more 
scatter, but below the freezing point most of their data lie 
above the lAPWS extrapolation. The measurements of Trinh 
and Ohsaka*'^ he below all other data; the deviation is about 
the same as their experimental uncertainty, which is 1%. 



down to about 250 K. Hacker's measurements show an inflec- 
tion point at about 268 K, below which the surface tension has 
a stronger temperattire dependence (Fig. 6). The data of Flo- 
riano and Angell are less accurate, but the authors also noted 
an inflection. The measurements of Trinh and Ohsaka show a 
systematic deviation from the other data, but the trend agrees 
with Hacker's data. Furthermore, a molecular dynamics sim- 
ulation by Lii and Wei*^ shows an inflection point as weU. 
LAPWS has recommended an equation to represent the sur- 



E. Speed of sound 

The speed of sound in supercooled water has been mea- 
sured in a broad frequency range, from 54 kHz to about 
20 GHz. In stable water at atmospheric pressure, no significant 
dependence on the frequency is found in this range, but the 
speed of sound in supercooled water shows a dispersion that 
increases with cooling. In this article we consider the thermo- 
dynamic (zero-frequency) limit of the speed of sound, which 
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is associated with the adiabatic compressibility.'^ Taschin et 
al.^'^ estimate that dispersion effects become noticeable at fre- 
quencies of 1 GHz and higher. We will therefore not con- 
sider the measurements in the 1-10 GHz range obtained with 
BriUouin light scattering, which show a speed-of-sound mini- 
mum between 250 K and 280 K; the temperature of the min- 
imum increases with increasing frequency."*^'^^ The lower- 
frequency (ultrasonic) measurements are shown in Fig. 7. Be- 
low 255 K, there appears to be dispersion; however, the speed 
of sound decreases with increasing frequency, which suggests 
an anomalous (negative) dispersion. While a negative disper- 
sion of the speed of sound very close to the vapor-liquid criti- 
cal point is, in principle, possible due to the divergence of the 
thermal conductivity,^" there is not yet any experimental in- 
dication for an anomaly of the thermal conductivity of super- 
cooled water.^' Recent simulations of Kumar and Stanley^^ 
show a minimum for this property. Debenedetti^ argues that 
negative dispersion can be ruled out because the BriUouin ex- 
periments all show positive dispersion. According to Taschin 
et al.f'^ the apparent negative dispersion is a result of sys- 
tematic errors in the data of Trinh and Apfel^' and Bacri and 
Rajaonarison.^^ The IAPWS-95 formulation agrees with the 
recent data of Taschin et al.^^ to within their uncertainty. 

There are no measurements of the speed of sound of super- 
cooled water at high pressure. In the stable region, Petitet et 
al.^^ performed measurements at 10 MHz up to 460 MPa, in- 
cluding the region close to the melting curve, down to 253 K. 
Petitet et al. have also published data down to 253 K at at- 
mospheric pressure, but these data were taken from Conde et 
al.^^ Since Conde et al. measured at 5 GHz, these data do not 
represent the zero-frequency limit of the speed of sound. 
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F. Spinodal 

In the supercooled region, IAPWS-95 predicts a re-entrant 
liquid spinodal, as shown in Fig. 8. The spinodal pressure be- 
comes positive at 233.6 K, which is a few degrees below the 
homogeneous nucleation limit. At about 195 K and 175 MPa, 
the spinodal curve crosses the homogeneous nucleation limit 
and enters the region where supercooled water can be experi- 
mentally observed. Up to about 290 MPa, the spinodal curve 
stays in the experimentally accessible range. Since a spinodal 
has not been observed there, the spinodal of IAPWS-95 con- 
tradicts experimental evidence. While a re-entrant spinodal is 
thermodynamically possible, Debenedetti has argued that the 
re-entrant spinodal scenario suggested by Speedy^^ is implau- 
sible for supercooled water. Debenedetti's argument may 
be summarized as follows. In the pressure-temperature plane 
(Fig. 8), a re-entrant spinodal must intersect the metastable 
continuation of the binodal, the vapor-pressure curve. At the 
intersection, liquid water is simultaneously in equilibrium 
with water vapor and unstable with respect to infinitesimal 
density fluctuations. Since these two conditions are mutu- 
ally exclusive, such an intersection cannot exist. On the other 
hand, at the vapor-liquid critical point the binodal and spin- 
odal do coincide, which is possible because the vapor and hq- 
uid phases are indistinguishable there. 



G. Liquid-liquid coexistence curve 

Both the existence of a second critical point and its location 
are still being debated in the Uterature. If the second critical 
point exists, there should be a liquid-liquid transition (LLT) 
curve - separating a hypothetical high-density hquid and low- 
density liquid - which ends at the critical point. At pressures 
below the critical pressure, water's response functions exhibit 
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Uquid-liquid critical point suggested by Bertrand and Anisimov.^ 



an extremum near the 'Widom line,' which is the extension 
of the LLT curve into the one-phase region and the locus of 
maximum fluctuations of the order parameter. 

While the location of the critical point obtained by differ- 
ent simulations varies greatly, different attempts to locate the 
LLT and the Widom line from experimental data have yielded 
approximately the same result. Kanno and Angell^^ fitted em- 
pirical power laws to their compressibility measurements and 
obtained singular temperatures located from 5 K to 12 K be- 
low the homogeneous nucleation temperature 7h (Fig. 9), sug- 
gesting a LLT that mimics the 7h curve but shifted to lower 
temperature. Mishima measured metastable melting curves of 
H2O ice IV" and D2O ices IV and V,^^ and found that they 
suddenly bent at temperatures of 4 K to 7 K below T^. Ac- 
cording to Mishima, this is indirect evidence for the location 
of the LLT, but a one-to-one correspondence between a break 
in the melting curve and the LLT has been questioned by Imre 
and Rzoska.^^ 

Mishima^^ approximated the LLT by a quadratic function 
of T with approximately the same shape as the Th curve, but 
allowing a shift to lower temperature. His final result, shown 
in Fig. 9, is close to Kanno and Angell's curve for pressures 
up to 100 MPa, albeit with a different curvature. 
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There are differences of up to 0.14% between the data 
sets, and it is not clear which is the best set. The isother- 
mal compressibility and thermal expansion coefficient have 
been measured in an extensive temperature and pressure 
range by Karmo and Angell.^^'^^ The isobaric heat capacity 
has been measured at atmospheric pressure by Angell and 
coworkers. ^'"^^ The speed of sound of supercooled D2O was 
reported by Conde et al.^^ at 5 GHz; at this frequency the 
speed of sound likely deviates from the zero-frequency limit, 
especially at lower temperatures. 



III. THERMODYNAMIC MODEL FOR SUPERCOOLED 
WATER 

In this section we further develop a scaling model for bulk 
thermodynamic properties supercooled water, which was ear- 
lier suggested by Fuentevilla and Anisimov' and, more re- 
cently, modified and elaborated by Bertrand and Anisimov.^ 



H. Heavy water 

Experimental data on the properties of supercooled 
heavy water (D2O) have been reviewed by Angell'^-'^ and 
Debenedetti.' The location of data for the density and its 
derivatives is shown in Fig. 10. At atmospheric pressure, the 
density has been measured by several investigators. •^^'^^■^^'^^ 



A. Scaling fields and thermodynamic properties 

Fluids belong to the universality class of Ising-Iike systems 
whose critical behavior is characterized by two independent 
scaling fields, a "strong" scaling field hi (ordering field) and 
a "weak" scaling field hi, and by a dependent scaling field 



7 



/z3 which asymptotically close to the critical point becomes a 
generalized homogeneous function of hi and /zi:^'^^^ 



,2-a-p 



(1) 



In this expression a ~ 0.1 10 and j3 ^ 0.326 are universal crit- 
ical exponents^^'^'* and with the superscripts ± referring 
to /j2 > and ^2 < 0, is a universal scaling function except 
for two system-dependent amplitudes. Associated with these 
scaling fields are two conjugate scaling densities, a strongly 
fluctuating scaling density (/)] (order parameter) and a weakly 
fluctuating scaling density (^2. such that 



d/j3 = 01 6h\ + <j>2 dh2 



with 



dhi 



h2 



dh2 



(2) 



(3) 



hi 



In addition one can define three susceptibilities, a "strong" 
susceptibility Xu a "weak" susceptibility X2, and a "cross" 
susceptibihty Xn- 
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In fluids and fluid mixtures one encounters a large variety of 
different types of critical phenomena.^^ The asymptotic ther- 
modynamic behavior near all kinds of critical points can be 
described in terms of Eq. (1). The differences arise from the 
actual relationships between the scaling fields and the physi- 
cal fields,'^ subject to the condition that at the critical point 



hi = h2 = h^ = 0. 



(6) 



In one-component fluids the relevant physical fields are the 
chemical potential n (Gibbs energy per mole), the tempera- 
ture T, and the pressure P. To satisfy condition (6) one defines 
= ^ - ^, AT = T - To, and AP = P - P^. In fliis article 
we adopt the usual convention that a subscript c refers to the 
value of the property at the critical point. There are two spe- 
cial models for critical behavior that deserve some attention. 
The first is the lattice gas in which the ordering field hi is 
asymptotically proportional to AjU and the weak scaling field 
proportional to AT?''~'^^ Hence, in the lattice gas 0i is pro- 
portional to Ap = p — po and 02 proportional to As = s — i^, 
where p is the molar density and s the entropy density. The 
lattice gas provides a model for the vapor-liquid critical point 
where the molar density yields the major contribution to the 
order parameter. In practice, the asymptotic critical behavior 
of a fluid near the vapor-Uquid critical point, including that of 
can be described by a slight modification of the lattice- 
gas model to account for some lack of vapor-liquid symmetry 
in real fluids. Another special model, introduced by Bertrand 
and Anisimov,^ is a "lattice liquid" in which the ordering field 



is asymptotically proportional to Ar and in which the weak 
scaling field is proportional to A^? Near the liquid-liquid 
critical point in weakly compressible supercooled water the 
entropy yields the major contribution to the order parameter 
and not the mass density, as first pointed out by Fuentevilla 
and Anisimov.^ Thus the thermodynamic properties near this 
liquid-liquid critical point can be described by a slight modi- 
fication of the lattice-Uquid model to account for some lack of 
symmetry in the order parameter.^ 

To implement a scaled thermodynamic representation it is 
convenient to make aU thermodynamic properties dimension- 
less in terms of the critical parameters 7^ and Pc or = p~^ : 



T = 



RTr 



PV, 



(7) 



where R is the ideal-gas constant. For the dimensionless phys- 
ical densities we define 



V = 



V 



s = 



R' 



Cp = 



R ' 



(8) 



where V is the molar volume, 5 the molar entropy, and Cp 
the isobaric molar heat capacity. The thermodynamic model 
of Bertrand and Anisimov was formulated in terms of P{jX,f) 
for which 



dP = V-^dfi + V-'SdT. 



(9) 



We have found it more convenient, and practically equivalent, 
to formulate this thermodynamic model in terms of fX(P,T) 
for which 



dfl=VdP-SdT. 



(10) 



Thus in this formulation, similar to that suggested earlier by 
Fuentevilla and Anisimov,^ we identify the order parameter 
with the entropy itself instead of the entropy density. In our 
model the scaling fields are related to the physical fields as 



hi =AT + a'KP, 
h2 = -AP + b'AT, 
h3 = AP - Aft + Afi'^ 



with 



AT = 



T-n 



AP = 



RTc 



Afl 



(11) 
(12) 

(13) 



(14) 



In Eq. (11) a' represents the slope —dt/dP of the phase- 
coexistence or Widom line at the critical point. In Eq. (12) b' is 
a so-called mixing coefficient which accounts for the fact that 
the critical phase transition in supercooled water is not com- 
pletely symmetric in terms of the entropy order parameter. In- 
troduction of mixing of this type is also known in the literature 
as the revised-scaUng approximation.^' Equation (1) only rep- 
resents the asymptotic behavior of the so-called singular criti- 
cal contributions to the thermodynamic properties. To obtain a 
complete representation of the thermodynamic properties we 
need to add a regular (i.e., analytic) background contribution. 
As has been common practice in developing scaled equations 
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of state in fluids near the vapor-liquid critical point, the 
regular background contribution is represented by a truncated 
Taylor-series expansion around the critical point: 

AA'^ £c™(Af)'"(AP)", with Coo = cio = coi = 0. 

. (15) 

The first two terms in the temperature expansion of AjX'^ de- 
pend on the choice of zero entropy and energy and do not 
appear in the expressions of any of the physically observable 
thermodynamic properties. Hence, these coefficients may be 
set to zero. Furthermore, the coefficient cqi = Vc — 1 = 0. 
Strictly speaking, critical fluctuations also yield an analytic 
contribution to h^.^^'^^ In this article we incorporate this con- 
tribution into the linear background contribution as has also 
been done often in the past. 

The current treatment of the background contribution 
differs from that in earlier publications.^'^ Previously, a 
temperature-dependent background was added to the critical 
part of each property. Because the background of each prop- 
erty was treated separately, the resulting property values were 
not mutually consistent. In this work, the background is added 
to the chemical potential, and the backgrounds in the derived 
properties follow, ensuring thermodynamic consistency. 

From the fundamental thermodynamic differential relation 
(10) it follows that 



5: 



dP/T 

dp. 
df 



= l-aVi+<|>2 + AA; 

^(j)i+b'(l)2-Afif. 



pi 



(16) 



(17) 



In this article we adopt the convention that a subscript P in- 
dicates a derivative with respect to P at constant f and a sub- 
script T a derivative with respect to T at constant P. Finally, 
the dimensionless isothermal compressibility kj, expansivity 
coefficient dy, and isobaric heat capacity Cp can be expressed 
in terms of the scaUng susceptibiUties Xi, Xi, and Xi2'- 



1 (dv\ 



(18) 
(19) 



[-a'xi+b'x2 + {'^-a'b')xn+^fifp] , 



B. Parametric equation of state 

It is not possible to write the scaled expression (1) for as 
an expUcit function of hi and h2. Such attempts always cause 
singular behavior of the thermodynamic potential in the one- 
phase region either at Ai = or at h2 = 0. This problem is 
solved by replacing the two independent scaling fields, hi and 
h2, with two parametric variables: a variable r which measures 



a "distance" from the critical point and an angular variable 
which measmes the location on a contour of constant r. A 
transformation most frequently adopted has the form: 



hi=ar^-"-Pe{l-e^), h2 = r{l~b^e^). 



(21) 



From Eqs. (1) and (3) it then follows that the order parameter 
<j)i must have the form:*'* 



<l>i=krl^M{e) 



(22) 



where M{9) is a universal analytic function of G. In princi- 
ple, this function can be calculated from the renormalization- 
group theory of critical phenomena.*^ In practice one adopts 
an analytic approximant for M(0), the simplest one being 
M(0) = 0:*6 



^i=kr^e. 



(23) 



Equations (21) and (23) define what is known as the "lin- 
ear model" parametric equation of state. In these equations a 
and k are two system-dependent amplitudes related to the two 
system-dependent amplitudes in Eq. (1), while b^ is a univer- 
sal constant which is often approximated by*' 



2-a-4j3 



(2-a-2/3)(l-2/3) 



1.361. 



(24) 



Equations (21) and (23) with the specific choice (24) for b^ 
is known as the "restricted" hnear model.^^ The resulting 
parametric equations for the various thermodynamic proper- 
ties can be found in many publications. ^^"*''**^'****'**^ In this 
paper we are using the "restricted" linear model to describe 
the supercooled-water anomaUes. The parametric equations 
needed for the analysis in this article are listed in Appendix A. 



IV. COMPARISON WITH EXPERIMENTAL DATA 

The scaling theory, formulated in the preceding section, 
represents the thermodynamic behavior asymptotically close 
to the critical point. Specifically, the liquid-liquid (LLT) curve 
that follows from Eq. (11) is a straight line, while the actual 
LLT curve should exhibit curvature (as shown in Fig. 9). In 
this section we investigate the thermodynamic properties of 
supercooled liquid water in a range of pressures and temper- 
atures, where the asymptotic theory appears to be adequate. 
Issues related to nonasymptotic features of the scaUng theory 
will be addressed in Sec. V. Hence, we restrict the asymp- 
totic theoretical model to pressures not exceeding 150 MPa, 
where the LLT can be reasonably approximated by a single 
straight line. The slope of the LLT was constrained to val- 
ues that are close to the slopes of the curves of Kanno and 
Angell^* and Mishima^' in the range of MPa to 150 MPa. 
Specifically, the value of a' in Eq. (11) was restricted to the 
range of 0.065 to 0.090. Because the position of the LLT is 
not precisely known, the critical point was allowed to deviate 
up to 3 K from Mishima's curve. It was found that the results 
of the model were rather insensitive to the critical pressure 
Pc, so Pc was constrained to the value of 27.5 MPa obtained 
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TABLE I. Parameter values for H2O in the model of Sec. IV 



Parameter 


Value^ 




Parameter 


Value 




Tc/K 


224.23 




Cll 


1.5363 X IQ- 


■ 1 


Pc/MPa 


27.5 




C12 


-6.4879 X IQ- 


■3 


Pc/(kg m"^) 


948.77 




C13 


7.7090 X 10^ 


3 


a 


0.22924 




C20 


-3.8888 X 10" 




k 


0.37704 




C21 


1.7347 X 10^ 


-1 


a' 


0.090 




C22 


^6.4157 X 10" 


-2 


C02 


7.1779 X 10" 


-2 


C23 


-6.9850 X 10" 


-3 


C03 


-4.093 6 X 10- 


-4 


C30 


6.9813X 10- 


1 


C04 


-1.0996 X 10" 


-3 


C3I 


-1.1459 X 10" 


-1 


C05 


2.9497 X 10" 


-4 


C32 


7.5006 X IQ- 


-2 



^ Final digits of parameter values are given to allow reproducing the values 
of properties with the model but do not have physical significance. 



TABLE IL Parameter values for D2O in the model of Sec. IV 



Parameter 


Value 




Parameter 


Value 




Tc/K 


232.65 




Cll 


1.2828 X 10" 


1 


Pc/MPa 


32.29 




fl2 


-1.6267 X 10" 


3 


Pc/(kg m-3) 


1055.74 




Cl3 


9.5552x 10" 


-3 


a 


0.22924 




C20 


-4.4118 X 10° 




k 


0.37704 




C2I 


3.0002 X 10- 


■1 


d 


0.078757 




C22 


-9.7204 X 10" 


■2 


C02 


6.9072 X 10- 


-2 


C23 


-1.4402 X 10- 


-2 


C03 


1.7651 X 10- 


-4 


C30 


8.4968X 10- 


-1 


C04 


-1.4458X 10" 


-3 


C3I 


-2.7188X 10- 


1 


C05 


4.3335 X 10" 


-4 


C32 


1.441 8 X 10" 


■1 



by Bertrand and Anisimov.^ It was also found that a nonzero 
mixing coefficient b' did not significantly improve the fit, so b' 
was set to zero. This means physically that the liquid-critical 
behavior in supercooled water exhibits little asymmetry in the 
order parameter and is indeed very close to lattice-liquid be- 
havior. 

Changes in the third decimal place of the values of the criti- 
cal exponents a and j3 result in small density changes that are 
of the order of 0.1%. However, some of the density measure- 
ments for water are more accurate than 0.1%; for example, 
the accuracy of the data of Hare and Sorensen'^ and Sotani 
et al}^ is 0.01%. Therefore, the values of the critical expo- 
nents must be given with at least four decimal places. We 
have adopted the values of Pelissetto and Vicari^^ and have set 
a = 0.1 100 and j3 = 0.3265. The value for the molar mass of 
H2O (18.015 268 g/mol) was taken from Wagner and PruB^^ 
and the ideal-gas constant R (8.3144621 J mol~^ K~^) was 
taken from Mohr et al.'^'^ The molar mass of D2O (20.027 508 
g/mol) was taken from an lAPWS guideline.^' 

The number of terms in the background Afl'^ [Eq. (15)] was 
increased step by step until the experimental data could be ac- 
curately represented. The final background contains fourteen 
free parameters. The reason for the many background terms of 
higher order in temperature and pressure is that the response 
functions are second derivatives of the thermodynamic poten- 
tial. To obtain, for example, a backgroimd term in the com- 
pressibility of second order in pressure, it is necessary to have 



a fourth-order pressure-dependent term in the potential. The 
terms in the backgrounds for each property are at most third 
order in temperature or in pressure. We want to emphasize 
that the observed anomalies are indeed due to the critical part 
of the equation of state since the nonlinear contributions to the 
regular part are needed only when the maximum pressure con- 
sidered is higher than about 100 MPa. Besides the background 
parameters, there are five additional parameters to be deter- 
mined: the critical temperature Tc and volume Vc, the linear- 
model amplitudes a and k, and the slope of the LLT Une a'. 
As noted, the values of and a' were constrained to a limited 
range. 

The model was fitted to heat-capacity data of Archer 
and Carter^^ and lAPWS-95, expansivity data of Hare and 
Sorensen,'** lAPWS-95 and Ter Minassian,^^ compressibil- 
ity data of Speedy and Angell,^^ Kanno and Angell^^ and 
Mishima,^' density data of Hare and Sorensen,'* Sotani et 
al.}'^ IAPWS-95 and Mishima,'^^ and speed-of-sound data 
of Taschin et al.^^ We have made small adjustments to the 
data of Mishima as described in Appendix C. For all quan- 
tities except the heat capacity, values from IAPWS-95 were 
only used at 0.1 MPa and above 273 K. For the heat capac- 
ity, IAPWS-95 values in the range of 273 K to 305 K up 
to 100 MPa were used. (We considered heat-capacity values 
calculated from IAPWS-95 more reliable than high-pressure 
heat-capacity data of Sirota et al?'') To reduce the time needed 
for optimization, not all data points were used in the fitting 
process; about 60 points were selected for each of the quan- 
tities heat capacity, expansivity, compressibility, and density. 
The selected data are given in the supplemental material. 

The model was optimized by minimizing the sum of 
squared residuals, where the residual is the difference be- 
tween model and experiment, divided by the experimental 
uncertainty.^^ For some data the uncertainty was not given 
and had to be estimated. The resulting optimized parameters 
are listed in Tables I and II. The value of a' for H2O is exactly 
0.09 because a' was restricted to the range of 0.065 to 0.090, 
and the optimum is located at the edge of this range. The fitted 
models are valid up to 300 K and from to 150 MPa. 

As can be seen in Fig. 11, the model represents the exper- 
imental density data well. The density jumps at low tempera- 
ture because the isobars cross the LLT curve there. In Fig. 12, 
the temperature of maximum density is plotted as a function 
of pressure, both for the model and for the IAPWS-95 formu- 
lation. At pressures higher than about 60 MPa, the values of 
IAPWS-95 deviate from the experimental data, while the cur- 
rent model agrees with the data (except at negative pressures, 
where the model was not fitted to any experimental data). The 
curve of the temperature of maximum density does not inter- 
sect the homogeneous nucleation curve at its break point at 
about 200 MPa and 180 K. Such an intersection was expected 
by Dougherty. 

The agreement of the compressibility data with values of 
the model is shown Fig. 13. An interesting feature is the inter- 
section of the isobars of 0.1 MPa and 10 MPa at about 250 K. 
The experimental data do not confirm or rule out such an in- 
tersection because of the scatter and the lack of data below 
245 K. However, the intersection implies that the pressure 
derivative of the compressibility, {dKj /dP)T, is positive at 
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Temperature (K) 



FIG. 11. Densities of H2O according to the model (curves). The sym- 
bols represent experimental data of Mishima,^' Sotani et al}^ and 
Hare and Sorensen.'** The symbols for Mishima's densities on dif- 
ferent isobars are altematingly open and fiUed to guide the eye. 




200 220 240 260 280 



Temperature (K) 

This work a Ter Minassian ef a/. (1981) 

This work x Henderson and Speedy (1987) 

IAPWS-95 □ Sotani et al. (2000) 

o Caldwell (1978) o Mishima (2010) 

FIG. 12. Temperature of maximum density of H2O as a function 
of pressure according to the model of Sec. IV (thick solid curve), 
the model of Sec. V (dashed) and IAPWS-95 (dotted). Tm marks 
the melting curve^-' and its extension to negative pressures;''-' Tn de- 
notes the homogeneous nucleation limit.^^ Symbols represent exper- 
imental data.'^'^^'^^'^^'^^ The temperatures of maximum density for 
Mishima's data^' were determined by locating the maxima of fits to 
his density data. 



low temperature and ordinary pressures. 

Figure 14 shows experimental data for the expansivity co- 
efficient and the values predicted by the model for five pres- 
sures. The model follows the experimental data, contrary to 
the IAPWS-95 formulation. At 240 K, where the difference 
between Hare and Sorensen's data of 1986 and 1987 is largest, 
the expansivity predicted by the model Ues between them. 

Heat-capacity data are compared with the model's predic- 
tions in Figs. 15 and 16. In Fig. 15, it is seen that the model 
follows the data of Archer and Carter,^"* whereas IAPWS-95 
follows the data of Angell et al.,^^ to which it was fitted. 
However, the curvature of the 0. 1 MPa isobar of the model is 
slightly higher than that suggested by the data of Archer and 
Carter. Murphy and Koop^^ proposed a heat-capacity curve 
with a broader peak than that of our model, but with about 
the same maximum value. The predicted isochoric heat capac- 
ity Cy diverges, as will be discussed in Sec. VIC. Figure 16 
shows the heat capacity as a function of pressure. There is a 
systematic difference between the data of Sirota et al}^ in the 
stable region and the values of IAPWS-95, and the data of 
Sirota et al. were not selected for the fit of the current model. 
At 250 K and pressures above about 50 MPa, the model pre- 
dicts a smaller pressure dependence of the heat capacity than 
IAPWS-95. The pressure dependence of the heat capacity is 
thermodynamically related to the expansivity coefficient, and 
we have seen that the expansivity coefficient of IAPWS-95 
does not agree with experimental data at low temperature and 
high pressure (see Fig. 14). Therefore, differences between the 
heat-capacity values of the current model and IAPWS-95 are 
to be expected. 

The speed of sound predicted by the model is shown in 
Fig. 7. The prediction agrees fairly well with the experimental 
data. At 240 K, the model shows a minimum in the speed of 
sound whereas IAPWS-95 predicts a monotonically decreas- 
ing speed of sound with decreasing temperature. The diver- 
gence of the speed of sound of the model at low temperature 
is related to a stability limit as will be discussed in Sec. VIC. 

For heavy water, the model was fitted to heat-capacity data 
of Angell et al?^ and Zabransky et al.,'^^ expansivity data 
of Kell^^ and Kanno and Angell, compressibility data of 
Kanno and Angell,^*^ density data of Kell,^' Zheleznyi,^^ and 
Rasmussen and MacKenzie,^^ and speed-of-sound data of 
Chen and Millero"^*^ and Marczak. It has been found that in 
the equation of state for H2O and D2O near the vapor-Uquid 
critical point deviations from corresponding states only ap- 
pear in the analytic background contributions, while the am- 
plitudes of the scaling fields are identical in accordance with 
corresponding states. '"^'^"^ Thus also for the second critical 
point we assigned to the linear-model amplitudes a and k the 
same values for D2O as for H2O. 

Figures 17-22 show the predictions of the model for heavy 
water compared with experimental data on density, com- 
pressibility, expansion coefficient, heat capacity, and speed 
of sound. The results are similar to those for ordinary water. 
The quality of the description is almost as good as for ordi- 
nary water, except for the speed of sound; the experimental 
data for speed of sound show noticeable (a few percent) de- 
viation from the model. The expansivity data of Zheleznyi,^^ 
shown in Fig. 20, deviate from the prediction of the model 
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FIG. 13. Isothermal compressibility of H2O according to the model 
(curves). For clarity, the curves are not shown for temperatures below 
the LLT line in the bottom graph. Symbols represent experimental 
data of Speedy and Angell,^^ Kanno and Angell,^^ and Mishima.^' 
Solid and open symbols with the same shape correspond to the same 
pressure. 



below about 255 K. For ordinary water, Hare and Sorensen'^ 
concluded that the data of Zheleznyi "show smaller densities 
for T < —25 °C than our bulk data and thus overestimate the 
anomaly in the expansivity." For heavy water, the density and 
expansivity values of Zheleznyi are also likely too low, since 
the densities measured by others^'^'^^'^^ are all higher. 

While the critical parts of the thermodynamic properties 
of H2O and D2O follow the law of corresponding states (the 
critical amplitudes a and k are the same) the corresponding 
regular parts do not follow this law. The critical compress- 
ibility factors Zc = Pc/{PcRTc) for H2O (Zc = 0.28) and D2O 
(Zc = 0.32) LLCP differ by about 13% while for the vapor- 
liquid critical point these factors differ by only \%}^^ The 
difference of 13% should be taken too seriously since the fit 
of the model is insensitive to the value of Pr at the LLCP. 
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FIG. 14. Expansivity coefficient of H2O according to the model 
(solid curves) and IAPWS-95 (dashed: within region of validity, dot- 
ted: extrapolations). Symbols represent experimental data of Ter Mi- 
nassian et al. and Hare and Sorensen. ■ 
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FIG. 15. Isobaric and isochoric heat capacity of H2O versus temper- 
ature at 0. 1 MPa according to the model (solid curves), IAPWS-95 
(dashed), and the prediction of Murphy and Koop (dotted). Sym- 
bols represent experimental data of Angell et al?^ and Archer and 
Carter.-''* The dash-dotted line is the estimated regular part of Cp. The 
predicted thermodynamic behavior of Cy is discussed in Sec. VIC. 
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FIG. 16. Isobaric heat capacity of H2O versus pressure according to 
the model (sohd curves) and IAPWS-95 (dashed). Symbols represent 
experimental data of Sirota et al?^ 
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FIG. 18. Temperature of maximum density of D2O as a function 
of pressure according to the model (thick soUd curve). Jm marks the 
melting curve^ and Th denotes the homogeneous nucleation Umit.^^ 
Symbols represent experimental data.^^'^^ 



a 
Q 




1040 



Zheleznyi (1969) 

Rasmussen and MacKenzie (1973) 
Emmet and Millero (1975) 
Kanno and Angell (1980) 
Hare and Sorensen (1986) 



220 240 260 280 
Temperature (K) 



300 



FIG. 17. Density of D2O according to the model (curves). The sym- 
bols represent experimental data.^^-^^"^^-^' 
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FIG. 19. Isothermal compressibility of D2O according to the 
model (curves). Symbols represent experimental data of Kanno and 
Angell.28 



V. SEMI-EMPIRICAL EXTENSION OF SCALING MODEL 

As mentioned in the preceding section, the asymptotic the- 
oretical model implies a linear LLT Une. An attempt to in- 
clude an additional term accounting for curvature of the LLT 
curve has been made by Fuentevilla and Anisimov^ by intro- 
ducing a pressure-dependent coefficient a' in Eq. (11). How- 
ever, such a procedure yields terms proportional to in some 
response functions which do not vanish far away from the crit- 
ical point, where the critical part should not play a role any- 
more. In principle, this problem can be solved by including 
crossover from singular critical behavior asymptotically close 



to the critical point to analytic behavior far away from the 
critical point, as has been done in the equation of state for 
H2O near the vapor-hquid critical point, and has also been 
suggested by Kiselev.^-^ However, it turns out that the theo- 
retical model can represent all experimental data for H2O up 
to the maximum available pressure of 400 MPa, if we sim- 
ply remove any constraints on the slope of the LLT and on 
the critical parameters in the unstable region which cannot be 
measured experimentally. With the addition of only two back- 
ground terms, we were then able to fit almost all experimental 
data for supercooled water with this semi-empirical extension 
of the theoretical model. 
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FIG. 20. Expansivity coefficient of D2O according to the model 
(curves). Symbols represent experimental data.^^-^^-^^ The data of 
Zheleznyi^^ (plusses) likely overestimate the anomaly in the expan- 
sivity (see the text). 
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FIG. 22. Speed of sound of D2O at 0.1 MPa according to the model 
(solid curve). Symbols represent experimental data of Conde et alr'^ 
and the dashed line shows the correlation of Marczak.''" 
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FIG. 21. Isobaric and isochoric heat capacity of D2O versus tem- 
perature at 0.1 MPa according to the model (solid curves). Symbols 
represent experimental Cp data of Angell et al?^ (squares: bulk wa- 
ter, circles: emulsion) and the dotted curve shows the correlation of 
Zabransky et al.'^^ for Cp. The predicted thermodynamic behavior of 
Cv is discussed in Sec. VIC. 



After fitting the model to a selection of experimental data 
at pressures up to 400 MPa, the parameters listed in Table III 
were obtained. Some parameters are significantly different 
from the earlier values. In paiticular, the critical pressure is 
about a factor of two higher than in the previous section, 
and the critical temperature is 10 K lower Nevertheless, the 
extended model does represent almost all available experi- 
mental thermodynamic data for supercooled water A com- 
parison between the density values predicted by the model 
is presented in Fig. 23. The model reproduces the data, ex- 
cept for Mishima's points between 160 MPa and 300 MPa 



TABLE III. Parameter values for the extended model of H2O 



Parameter 


Value 




Parameter 


Value 




TJK 


213.89 




Cl2 


-4.6569 X 10- 


i 


Pc/MPa 


56.989 




Cl3 


2.3627 X 10^ 


3 


Pc/(kgm-3) 


949.87 




Ci4 


-2.8697 X 10^ 


4 


a 


0.11624 




C20 


-3.6144 X 10" 




k 


0.43280 




C21 


-1.5009 X 10" 


2 


a' 


0.10898 




C22 


-2.4609 X 10" 


2 


C02 


4.0793 X 10" 


-2 


C23 


9.8679 X 10" 


4 


C03 


-6.7912 X 10" 


-4 


C30 


5.4267 X 10" 


1 


C04 


-7.5669 X 10" 


-6 


C3I 


1.0620X 10" 


1 


C05 


1.0922 X 10" 


-5 


C32 


1.2759X 10" 


2 


C'll 


1.9547 X 10" 


-1 


C41 


-7.9970 X 10" 


2 



below 230 K. The temperature of maximum density is plot- 
ted in Fig. 12. The calculated behavior of the extended model 
is similar to that of the asymptotic model in the previous sec- 
tion in the range from MPa to 120 MPa. At higher pressures, 
the extended model predicts higher temperatures of maximum 
density than the previous model, but the results are still within 
the experimental uncertainty. 

Compressibility data are compared with values of the ex- 
tended model in Fig. 24. As in the previous model, the isobars 
of 0.1 MPa and 10 MPa intersect, but the intersection is lo- 
cated at a lower temperature. The data of Speedy and Angell^^ 
and Kanno and Angell^^ are well represented. The extended 
model also reproduces most of Mishima's data, which have a 
lower accuracy than the data of Angell and coworkers. Fig- 
ure 25 shows the expansivity coefficient predicted by the ex- 
tended model, which agrees with the data of Ter Minassian et 
al.^^ Below 250 K, the model agrees better with the data of 
Hare and Sorensen of 1986 than with their data of 1987. 

The calculated heat capacity at 0.1 MPa is compared with 
experimental data in Fig. 26. For the extended model, the 
maximum of the heat capacity is lower than for the asymptotic 
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FIG. 23. Densities of H2O according to the extended model (curves). 
The symbols represent experimental data of Mishima,^' Sotani et 
al}^ and Hare and Sorensen.'** The symbols for Mishima's densities 
on different isobars are alternatingly open and filled to guide the eye. 
The vertical Unes through Mishima's points are uncertainties given 
by Mishima. 



model, and the curvature of the data of Archer and Carter^'' 
is better represented. The predicted heat capacity is shown 
as a function of pressure in Fig. 27. There are large differ- 
ences between the results of the model and those of lAPWS- 
95; the model predicts a minimum in the 250 K heat-capacity 
isotherm at about 240 MPa. A minimum in the heat capacity 
at this location was also predicted by Ter Minassian et al?^ 
based on their measurements of the expansivity coefficient. 

The speed of sound predicted by the model, shown in Fig. 7, 
agrees fairly well with the experimental data of Taschin et 



VI. PHYSICAL INTERPRETATION OF THE MODEL 

A. Two states in supercooled water 

Two features make the second critical point in water phe- 
nomenologicaUy different from the well-known vapor-liquid 
critical point. The negative slope of the liquid-liquid phase 
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FIG. 24. Isothermal compressibility of H2O according to the ex- 
tended model (curves). Symbols represent experimental data of 
Speedy and Angell,^^ Kanno and Angell,^^ and Mishima.^' SoUd 
and open symbols with the same shape correspond to the same pres- 
sure. The regular backgrounds for the compressibility are shown by 
dash-dotted lines. The negative values of the regular parts are due 
to incorporating the fluctuation-induced critical background in Xl as 
noted in Appendix A. 
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FIG. 25. Expansivity coefficient of H2O to the extended model 
(curves). Symbols represent experimental data of Ter Minassian et 
al.^^ and Hare and Sorensen. '^'^^ The dash-dotted Une is the esti- 
mated regular part at 0. 1 MPa. 
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FIG. 26. Isobaric and isochoric heat capacity of H2O versus tem- 
perature at 0.1 MPa according to the extended model (solid curves), 
IAPWS-95 (dashed), and the prediction of Murphy and Koop (dot- 
ted). Symbols represent experimental data of Angell et al?^ and 
Archer and Carter. -^"^ The dash-dotted line is the estimated regular 
part of Cp. The predicted thermodynamic behavior of Cy is discussed 
in Sec. VIC. 
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FIG. 27. Isobaric heat capacity of H2O versus pressure according to 
the extended model (solid curves) and IAPWS-95 (dashed). Symbols 
represent experimental data of Sirota et alP 



transition line in the P-T plane means that high-density liq- 
uid water is the phase with larger entropy. The relatively large 
value of this slope at the liquid-liquid critical point (about 25 
times greater than for the vapor-liquid transition at the crit- 
ical point) indicates the significance of the entropy change 
relative to the density change, and, correspondingly, the im- 
portance of the entropy fluctuations. These features suggest 
that liquid-liquid phase separation in water is mostly driven 
by entropy rather than by energy, thus supporting the "lattice- 
liquid" choice for the scaling fields given by Eqs. (11)-(13) 

with y = 0. 



As the first step to understand a relation between water's 
polyamorphism and the behavior of cold aqueous solutions, 
Bertrand and Anisimov^ have introduced a mean-field "two- 
state" model that clarifies the nature of the phase separation 
in a polyamorphic single-component liquid. Two-fluid-states 
models trace their lineage back to a 19th century paper by 
Rontgen."''* Relatively recently, Ponyatovskii et al}^^ and, 
more quantitatively, Moynihan^"^ have described the emer- 
gence of a LLCP in supercooled water as resulting from the 
effects of nonideality in a mixture of two "components," with 
their fraction being controlled by thermodynamic equilibrium. 
However, while Moynihan assumed a "regular-solution" type 
of nonideality, which implies an energy-driven phase separa- 
tion, such as the vapor-liquid transition or the conventional 
liquid-liquid transition in binary solutions, we believe that a 
near "athermal-solution" type of nonideality is mainly respon- 
sible for the entropy-driven phase separation in metastable 
water near the LLCP. 

It is assumed that liquid water is a "mixture" of two states, 
A and B, of the same molecular species. For instance, these 
two states could represent two different aiTangements of the 
hydrogen-bond network in water and correspond to the low- 
density and high-density states of water The fraction of water 
molecules, involved in either structure, denoted for state A 
and 1 — for state B, is controlled by thermodynamic equi- 
librium between these two structures. Unlike a binary fluid, 
the fraction is not an independent variable, but is deter- 
mined as a function of pressure P and temperature T from the 
condition of thermodynamic equilibrium. The simplest "ather- 
mal solution model"'"^ predicts a symmetric liquid-liquid 
phase separation for any temperature, with the critical frac- 
tion 0c = 1 /2, if the interaction parameter, which controls the 
excess (non-ideal) entropy of mixing is higher than its criti- 
cal value. However, unlike a purely athermal non-ideal binary 
fluid, the entropy-driven phase separation in a polyamorphic 
single-component liquid will not be present at any tempera- 
ture. On the contrary, the critical temperature is to be spec- 
ified through the temperature dependence of the equilibrium 
constant K by 

ln^ = A(l-l), (25) 

where A is the heat (enthalpy change) of "reaction" between 
A and B. A finite slope of liquid-liquid coexistence in the P- 
T plane can be incorporated into the two-state lattice liquid 
model if one assumes that the Gibbs energy change of the "re- 
action" also depends on pressure. 

A conceptually similar two-state model suggested by Hruby 
and Holten'"^ has explained the inflection point in the surface 
tension (Fig. 6) by a rapid increase of the fraction of water 
molecules in the low-density state as water is cooled down. 

B. Liquid-liquid phase separation 

A consequence of the second critical point is a phase sepa- 
ration into a high-density and a low-density liquid at tempera- 
tures below the critical temperature. From the "lattice-liquid" 
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FIG. 28. Temperature-density diagram. Ci and C2 indicate ttie first 
and second critical point. The line left of Ci is the saturated vapor 
density according to IAPWS-95. Down to 250 K, the line marked 
by 'Liquid' is the saturated liquid density computed with lAPWS- 
95; below 250 K, the line represents the liquid density at 0.1 MPa 
predicted by our model. The dashed lines show the densities of the 
two liquid phases in equilibrium on the LLT line. 
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FIG. 30. Different linearizations of the LLT curve. The thick solid 
curve is the LLT curve suggested by Mishima.^' The lines marked 
a and b are the linearized LLT lines of our models restricted to 
150 MPa and 400 MPa, respectively; the line marked c is the LLT line 
of Moynihan.'"^ Symbols mark the liquid-liquid critical points; the 
solid diamond is the critical point of Mishima,^' the solid square is 
that of Bertrand and Anisimov,-' the open circle is that of our asymp- 
totic model, the crossed circle is that of our extended model, and the 
cross is that of Moynihan.'^^ The dotted line is the phase-transition 
line between the two amorphous phases LDA and HDA as estimated 
by Whalley et al.^^° 
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FIG. 29. Temperature-density diagram of the two liquid phases in 
equilibrium on the LLT line. Solid line: asymptotic model, dashed 
line: extended model, dotted line: prediction of the two-state regular- 
solution model of Moynihan.'"^ Solid dots mark the second criti- 
cal points. Open circles are the densities of the low-density amor- 
phous (LDA) and high-density amorphous (HDA) phases of water 
estimated from the experiments of Mishima.'"^ 



two-state model with a steep slope of the liquid-liquid tran- 
sition (LLT) line in the P-T phase diagram it follows that 
the phase separation is mostly driven by entropy rather than 
density. The LLT line is not exactly vertical in the phase dia- 
gram due to a small density difference between the two liquid 
phases, as pointed out by Bertrand and Anisimov.^ In Fig. 28, 
a temperature-density phase diagram is shown including both 
the vapor-liquid and liquid-liquid critical points. The region 



below the second critical point is shown in more detail in 
Fig. 29, where predictions by several models of the densities 
of the two liquid phases are compared. 

In the second-critical-point scenario it is assumed that the 
LLT line is connected to the phase-transition line between 
two amorphous phases found below about 130 K.^'^ Since 
the amorphous phase-transition line is nearly horizontal in the 
P-T diagram, as shown in Fig. 30, it follows that the LLT 
must be strongly curved to account for such a connection. 
In our model, we have linearized the LLT, so the predicted 
equilibrium densities at 135 K differ from the experimental 
densities of the amorphous phases, as seen in Fig. 29. The 
two-state model of Moynihan"^^ was fitted to the experimen- 
tal amorphous-phase densities, as shown in Fig. 29. In Moyni- 
han's model, the LLT is also a straight line, but its slope, 
dP/dr, is less than half of the slope of our LLT lines (Fig. 30). 
Comparing our models for the 0-150 MPa and 0^00 MPa 
ranges, we see that a lower LLT slope is needed for the larger 
pressure range. In conclusion, with a linearized LLT, the larger 
the desired temperature and pressure range of the model, the 
smaller dP/dT slope of the linearized LLT line is required. 

Based on the fit of the model to the experimental data, the 
liquid-liquid critical point can be located in a certain range of 
temperatures and pressures, depending on the pressure limit 
for experimental data involved in the fit. Without any prelim- 
inary idea on the location of the LLT curve, the position of 
the liquid-liquid critical point is uncertain. In principle, we 
cannot exclude even a negative value for the critical pressure. 
Our assumption on the LLT curve is based on the predic- 
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FIG. 31. Reduced sum of squared residuals as a function of the location of the liquid-liquid critical point for H2O. (a) Asymptotic model, 
fitted to experimental data up to 150 MPa; (b) Extended model, fitted to experimental data up to 400 MPa. 



tions of different authors^^'^^'^^'^'^ which give approximately 
the same result. We have linearized this curve for a certain 
range of pressures. If the data are restricted up to 100 MPa, 
the critical pressure is optimized to 28 MPa and the criti- 
cal temperature is found to be 227 K.^ If we include more 
high-pressure data and linearize the LLT curve in the pressure 
range up to 150 MPa, the optimal value for the critical pres- 
sure becomes less certain and moves up. This tendency con- 
tinues when we include high-pressure data and linearize the 
LLT cm-ve up to 400 MPa. As illusti-ated in Fig. 31(a), best 
fits for the 150 MPa model are obtained with critical-pressure 
values between about 25 MPa and 70 MPa. For the extended, 
400 MPa, model, the optimum critical-pressure values move 
to a higher range from 40 MPa to 85 MPa [Fig. 31(b)]. 



C. Absolute stability ilmit of the liquid state 
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For the supercooled liquid state to be thermally and me- 
chanically stable, both the isochoric heat capacity Cv and the 
isothermal compressibility Kj must be positive. The scaling 
model predicts regions in the phase diagram where Cy or Kt 
is negative, as shown in Fig. 32. As noted by Bertrand and 
Anisimov,^ the stability locus is located close to the Widom 
line, but its exact location depends on the parameters of the 
model. Our model predicts that the thermal stability condi- 
tion is violated before the mechanical stability condition as 
water is cooled down, contrary to the violation order shown 
by a cubic equation of state that describes the vapor-liquid 
transition. ^'^'''^ At pressures slightly above the critical pres- 
sure, the mechanical stability limit coincides with the LLT 
line; that is, one of the liquid phases (the low-density liquid) 
is unstable. We note that a spinodal for the LLT is not defined 
in the scaled equation of state.' 



FIG. 32. Absolute stability limit of the liquid state, predicted by the 
model. Solid line: limit of mechanical stability, where the isothermal 
compressibility Kj is zero. Long dashed line: limit of thermal stabil- 
ity, where the isochoric heat capacity Cy is zero. Th marks the exper- 
imental homogeneous nucleation limit, and C2 indicates the second 
critical point. The dotted line, at which hi = [Eq. (11)], indicates 
the liquid-liquid transition line and its extension below the critical 
pressure, the Widom line. The inset shows the stability limits in the 
vicinity of the critical point. 



The isochoric heat capacity Cy is related to the isobaric heat 
capacity by' 

Cv^Cp-^^. (26) 
PKt 
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FIG. 33. Contribution of the susceptibilities, Xi, Xi^ and and the 
regular part to the compressibiUty, as given by Eq. (18), at 10 MPa. 
The thick line represents the total compressibiUty multiplied by the 
volume, kjV . 

The critical part of expressed through the scaling suscepti- 
bilities is derived in Appendix B. As the compressibility ffr 
becomes smaller with decreasing temperature, Cy first be- 
comes negative and then diverges to negative infinity as Kj ap- 
proaches zero. At temperatures below the mechanical stability 
limit, where Kj is negative, Cy is larger than Cp (Fig. 15). An- 
other property that diverges is the speed of sound c, which is 
related to Cp, Cy and fCr according to^ 

C={pKTCylCpY^I^. {21) 

The speed of sound diverges to positive infinity (Fig. 7) 
when Cy approaches zero (at the thermal stability limit), 
which occOTS at a slightly higher temperature than the zero- 
compressibility limit. In the unstable region, the speed of 
sound is imaginary. 

Kiselev^ has computed a kinetic spinodal of supercooled 
water and found that this stability limit was located a few 
kelvin below the homogeneous nucleation limit. Kiselev calls 
the area in the phase diagram below the stability limit a "non- 
thermodynamic habitat," not because liquid water is neces- 
sarily mechanically unstable there, but because the hfetime 
of liquid water is smaller than the time needed for equili- 
bration. Using molecular simulations, Moore and Molinero"*' 
have also found a stability limit just below the homogeneous 
nucleation temperature. They show that below this Umit, the 
time for crystallization is shorter than the relaxation time of 
the liquid, so liquid water cannot not be equilibrated. 

An important question arises: What terms in our model are 
responsible for the instability of the liquid state? Figure 33 
shows the contributions of the individual susceptibilities to the 
compressibility, as given by Eq. (18). The strong and weak 
susceptibiUties Xi and Xi are positive in the entire tempera- 
ture range, whereas the cross susceptibility Xi2 changes sign 
on the Widom line. The sum of the contributions of the three 
susceptibiUties is positive in the entire range; it is the contri- 
bution of the regular part, which is negative, that results in a 
negative compressibility below a certain temperature. 
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FIG. 34. Experimental densities of D2O (open and filled cir- 
cles, both in hydrophiUc confinement),''^'"^ and H2O (open 
squares: hydrophobic confinement, filled squares: hydrophilic 
confinement),"*^'"^ together with bulk densities '^'*'^ (plus and cross 
symbols). The predictions of the asymptotic scaling models for bulk 
D2O and H2O are shown by solid curves, and the dashed curves rep- 
resent the extrapolation in the unstable region. 



D. Confined water 

New insights into the thermodynamics of supercooled wa- 
ter have been provided by measurements of water confined in 
mesoporous structures. The siUca material MCM-41, which is 
comprised of a "honeycomb" of long cylindrical nanopores, 
has received significant attention as a confining medium in 
part because the diameter of the pores can be precisely con- 
trolled. As the pore diameter is decreased, the freezing point 
of water is reduced. Spontaneous crystallization is completely 
suppressed below a critical pore diameter, which allows for 
measurements to be made below the homogeneous-nucleation 
temperature. Examples of such measurements for hydrated 
MCM-41 are presented for the density in Fig. 34, for the ex- 
pansion coefficient in Fig. 35, and for the heat capacity in 
Fig. 36, to be compared with the prediction for bulk water. 

Many of the results of experiments on confined water can 
be interpreted in the framework of the second critical point 
hypothesis. Remarkably, the heat capacities exhibit maxima 
upon crossing the Widom line predicted for bulk water. How- 
ever, the height of the maximum and the shape of the heat 
capacity in porous media differ significantly from those pre- 
dicted in bulk water (Fig. 36). It has recently been reported 
that the density of confined water shows significant hysteresis 
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FIG. 35. Thermal expansion coefficient in confined water as mea- 
sured by MaUamace (crosses).^^" For comparison, the expansion co- 
efficient of bulk water is also shown (circles). The curve is the pre- 
diction of the asymptotic model. 
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FIG. 36. Heat capacity of H2O and D2O in silica nanopores as mea- 
sured by Nagoe et al. Different symbols correspond to different 
pore diameters. The curves are the predictions of the asymptotic 
model. 



with temperature variation at elevated pressures. This can 
be interpreted in terms of the crossing of the first-order Uquid- 
liquid coexistence curve. The exact connection between the 
bulk thermodynamics of supercooled water and the results 
found in confinement is an active area of research. Surface ef- 
fects, both in the form of finite-size effects and surface-water 



interactions, are thought to be contributing factors to the ob- 
served differences. In computer simulations of water confined 
in MCM-41, investigators have found that water found in the 
inmiediate vicinity of the surface behaves differently from wa- 
ter in the center of the pore.'^'^ This distinction raises the pos- 
sibility that wetting-type phenomena may be responsible for 
some of the observed behavior. That the surface interactions 
play a significant role can be seen from the fact that the mea- 
sured density for hydrophobic and hydrophilic confinement 
are qualitatively different.''^ However, even in similar hy- 
drophiUc environments, the measurements for D2O show con- 
siderable difference (Fig. 34). The difference could be caused 
by the uncertainty in the background of the neutron diffrac- 
tion intensity, which must be subtracted from the measured 
intensity in neutron- scattering experiments. ^'•^'^'^ In addition, 
the relative densities that are obtained must be scaled to ob- 
tain absolute densities, which introduces an additional uncer- 
tainty. The shape of the response functions, the expansivity 
coefficient and the heat capacity, as shown in Figs. 35 and 36, 
is more certain. The rounded anomalies of the heat capacity 
around the Widom line in the confined mW water model are 
demonstrated by recent simulations of Xu and Molinero.^^^ 
An initial qualitative analysis of the finite-size effects based 
on the LLCP hypothesis has been made by Bertrand and 
Anisimov.^ This work indicates that finite-size effects can ac- 
count for some of the differences between confined and bulk 
results. 



VII. DISCUSSION 

We have evaluated the IAPWS-95 formulation for ordinary 
water when it is extrapolated into the supercooled region. At 
0.1 MPa, it was found that IAPWS-95 reproduces the exper- 
imental data of all properties except the isobaric heat capac- 
ity. With increasing pressure, however, the deviations of the 
IAPWS-95 values from the data become larger and larger. A 
property that is particularly poorly described is the expan- 
sion coefficient, for which the IAPWS-95 extrapolation has 
a wrong sign in a large region of the phase diagram. Even 
in the stable region close to the melting line, where IAPWS- 
95 should be valid, the expansion coefficient predicted by 
IAPWS-95 is a factor of two lower than the experimental 
value. The lAPWS equation for the surface tension performs 
well when extrapolated down to about 266 K; below that, the 
experimentally observed increase of the temperature depen- 
dence is not reproduced and the difference between the equa- 
tion and the data increases with decreasing temperature. 

We have demonstrated that a theoretical model based on 
the assumption of a liquid-liquid critical point in supercooled 
water can represent the thermodynamic properties of both su- 
percooled H2O and D2O to pressures of 150 MPa. Moreover, 
by allowing the slope of the liquid-liquid transition (LLT) line 
and the critical pressure to be freely adjustable parameters the 
model can represent almost all available thermodynamic prop- 
erty data for supercooled water. Nevertheless, there are stiU a 
number of issues that need to be considered. 

A principal issue is that the existence of a liquid-Uquid crit- 
ical point is not the only possible explanation for the anoma- 
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lous behavior of the thermodynamic properties of supercooled 
water. Scenarios for a singularity-free or critical-point-free in- 
terpretation have also been proposed. '^'^"^^^ Recently, Stokely 
et al. have shown, for a water-like lattice model, how all 
such scenarios can be described by varying two quantities, 
the strength of the hydrogen bonds and the cooperativity of 
the hydrogen bonds. An intriguing possibility of the existence 
of multiple critical points in supercooled water, as predicted 
by some simulations,'^**''^' is discussed by Brovchenko and 
Oleinikova.'^*^ Another possibility is that response functions 
like the compressibility do not diverge at a single temper- 
ature corresponding to a critical temperature but at a range 
of pressure-dependent temperatures Ts (P) corresponding to a 
crystallization spinodal, the absolute limit of stability of the 
liquid phase. Most recently, the discussion on the nature of 
the anomalies observed in supercooled water received an addi- 
tional impetus after Limmer and Chandler reported new sim- 
ulation results' for two atomistic models of water, mW^^^ and 
mST2.^'^ They did not find two liquid states in the super- 
cooled region and excluded the possibility of a critical point 
for the models studied. However, the most recent simulations 
by Sciortino et al. of the original ST2 model confirm the 
existence of a liquid-liquid critical point for that particular 
model, with two distinct liquid states. They found no evidence 
of crystallization during their simulation time, while there was 
enough time for the liquid to equiUbrate. It would be impor- 
tant to compare the anomalies obtained for various models 
with those exhibited by real water. 

Recent simulations of Moore and MoUnero^'^ indicate that 
for simulation cells larger than the critical ice nucleus, spon- 
taneous crystallization occurs before liquid-liquid separation 
can eqiuhbrate. An approach to consider a couphng between 
spontaneous crystallization near the absolute stability limit of 
the liquid state and liquid-liquid separation is the theory of 
weak crystallization pioneered by Brazovskii,'^"* further de- 
veloped by Brazovskii et al. and reviewed by Kats et al. 
According to this theory, the fluctuations of the translational 
(short-wavelength) order parameter y/ renormalize the mean- 
field distance A(, — {T — Tq)/T between the temperature T and 
the mean-field absolute stabiUty limit 7b of the liquid phase: 

A = Ao + ffA-i/2, (28) 

where A = {T — Ts)/T where Ts is the fluctuation-affected 
stability-limit temperature, and JC is a molecular parame- 
ter, similar to the Ginzburg number that defines the valid- 
ity of the mean-field approximation in the theory of critical 
phenomena. **^''^^ The theory of weak crystallization requires 
fcA^'/^ <C Ao and Aq <C 1. There are two non-trivial con- 
sequences of the effects of translational-order fluctuations: 
(1) The heat capacity would contain a fluctuation correction 
5Cp ~ A~^/^; (2) The renormalized gap A may be positive 
even at negative Aq, meaning that translational-order fluctua- 
tions stabilize the liquid phase. 

The theory of weak crystallization was used to de- 
scribe a coupling between one-dimensional density mod- 
ulation and orientational fluctuations in liquid crystals; it 
was supported by accurate light-scattering and heat-capacity 
measurements. ^^^"^^' 



Similarly, we could consider a coupling between the 
translational-order fluctuations and the critical fluctuations of 
the liquid-liquid order parameter . A relationship between 
the structural transformation in Uquid water and its crystal- 
lization rate is suggested by recent simulations by Moore 
and Molinero."^ The lowest-order term associated with such 
a coupling in the free-energy Landau expansion would be 
~ Such a coupling would be possible even if the virtual 
critical point is below the stability limit of the liquid phase. 

Formation of ice is not "weak crystalUzation." This is why 
application of weak crystallization theory to supercooled wa- 
ter is questionable. However, it may not be hopeless. The gap 
A between the melting temperature and the temperatiu'e of 
spontaneous crystallization 7h is of the order of 0.1. Hence, a 
further investigation of the possibility of such coupling would 
be worthwhile. The actual question is how the experimen- 
tal observations can be explained by a theory that accounts 
for both crystalUzation and liquid-liquid separation. The fi- 
nal conclusion concerning the existence of the liquid-liquid 
critical point in water should be based on the ability to quan- 
titatively describe the experimental data. 

Assuming the existence of the LLCP in supercooled wa- 
ter, we confirm the finding of Fuentevilla and Anisimov^ and 
Bertrand and Anisimov^ that the critical point is located at 
much lower pressure than predicted by most simulations, def- 
initely below 100 MPa, while the precise value of the critical 
pressure is uncertain. However, regardless of all open ques- 
tions, we have shown that a critical-point parametric equation 
of state describes the available thermodynamic data for super- 
cooled water within experimental accuracy, thus establishing 
a benchmark for any further developments in this area. 
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Appendix A: Linear-model parametric equations 

The formulas for the linear model in this appendix are taken 
from Behnejad et a/.^' [In Behnejad et al.,^^ the formula cor- 
responding to our Eq. (A15) erroneously does not contain the 
factor Lo-] 
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The r and 9 coordinates are related to h\ and h2 by 



hi=ar^+y6{\-e'^) 
h2 = r 

with7=2-a-2j3 and 



h2 = r{\-b^6^), 



r(i-2j3)- 

The scaling densities are given by 

<j>2 = akr^-"s{d), 

with 

S2 = ia-l){y-2li)b\ 
The scaling susceptibilities are given by 

k 



(Al) 
(A2) 



(A3) 



(A4) 
(A5) 

(A6) 
(A7) 
(A8) 
(A9) 



Xi 



(AlO) 
(All) 



Xi2= krP-'qn{6), 
Xi=akr-"q2{e), (A12) 

with 

qi{6)= {l-b^e^ + 2l5b^e^)/qo{e), (A13) 
^12(0)= 0[-r+(r-2/3)02]/^o(0), (A14) 
^2(0)= [(l-a)(l -302)^(0) 

- (j3 + r)2s2L„eHi - 0')] /qo{e), (ai5) 

^o(0) = {l-3d^){l-b^d^) + 2b\l5 + Y)d\l - 02). 

(A16) 

The weakly divergent scahng susceptibility Xi also contains 
a negative fluctuation-induced analytical background, which, 
as suggested by Bertrand and Anisimov,-' was incorporated in 
the regular part of the chemical potential in our approach. 

Appendix B: Isochoric heat capacity in scaling models 



Equations (11)-(13) can be generalized as 

hi = cAf + aAP, 
h2 = -dAP + bAt, 
h3 = AP - Aft + Afl\ 



(Bl) 
(B2) 
(B3) 



The primes of the parameters a' and /?' have been omitted for 
brevity. With a = b = and c = d = I we recover the pure 
lattice-hquid model, whereas with a = b = 1 and c = d = 
we obtain the pure lattice-gas model. The case c = for the 



lattice-gas model corresponds to hi = AP = Afi, as discussed 
in Ref. 3. The expression for the isochoric heat capacity Cv, 



Cv =Cp-Vf^. 

Kt 



(B4) 



leads to the critical part of Cy, 



Cv ZlZ2-Z?2 

f a^Xi-'^adXn+d^Xi' 



(B5) 



with A = {ab + cd)^. For the pure lattice liquid, we obtain 



t Xi' 



and for the pure lattice gas. 



(B6) 



(B7) 



as obtained by Behnejad et al.^^ In our model, the coefficient 
a is small but finite and Cy is given by 



(B8) 



Cv ^ XxXi-Xn 
f a^Xi-'^Xn+X2 

Appendix C: Adjustment of experimental data 



The density data of Mishima show systematic differences 
of up to 5 kg/m-' with the densities reported by Asada et 
al.,^^ with lAPWS-95 (in the region where it is reliable), and 
with the equation of state of Saul and Wagner.''*'^ The dif- 
ferences appear to be mostly pressure dependent and only 
slightly temperature dependent; see Fig. 37. In the range from 
273 K to 373 K and MPa to 380 MPa, where the three other 
data sets (Asada et al., lAPWS-95, Saul and Wagner) over- 
lap and can be considered rehable, their maximum mutual 
difference is 0.8 kg/m'. Therefore, we decided that a correc- 
tion of Mishima's data is justified. A pressure-dependent den- 
sity correction was determined by fitting a quadratic function 
of the pressure to the difference of Mishima's densities and 
the correlation'^' of Asada et al., between 245 K and 274 K 
(Fig. 37). This density correction was then subtracted from all 
densities measured by Mishima. 

Because of an incorrect graph, the heat-capacity data shown 
in the article of Bertolini et al.^^^ are about 4% too low. The 
discrepancy was noted by Duan et al.,^'^^ but they did not ap- 
ply a correction. Tombari et al.^^ and Mishima^' corrected 
the data without comment. A quantitative correction can be 
obtained because the graph of Bertolini et al. also contains a 
plot of an equation of Leyendekkers and Hunter,'"'^ and the 
data in the graph can be scaled until the curve matches the 
equation. The corrected heat capacity C|,°" may be calculated 
from the original C"^ by C^p°" = 13.74 0.863 C^^ with Cp 
in Jmol~^K~'. This correction is slightly smaUer than that of 
Tombari et al. and of Mishima. 

'p. G. Debenedetti, J. Phys.: Condens. Matter 15, R1669 (2003). 
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FIG. 37. Difference of Mishima's^^ density and the correlation of 

Asada et al}^ The curve is a quadratic fit, Ap = -4.83 + 0.03348/"- 
3.87 X IQ^^P^, with Ap in kg/m^ and pressure P in MPa. 
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